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We study the superconducting transition in the attractive Hubbard model in two dimensions using 
the self-consistent T-matrix approximation. We demonstrate that for large system sizes, this ap- 
proximate method produces XY critical scaling in the correlation length and pair susceptibility. For 
the parameters we investigate, the critical regime is quite large, extending beyond 5 times Tktb- We 
discuss the role of vortex-pair unbinding in the normal state in the context of pseudogap behaviour 
and the relevance to the physics of underdoped cuprates. 
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Compared to conventional superconductors, the ratio 
of the superconducting pair wave- function size to the av- 
erage inter-carrier spacing is many orders of magnitude 
lower in the cuprates. Consequently, the pairing in the 
cuprates is much less mean-field like and quite some time 
ago it was argued [Q that, in the cuprates, strong pairing 
correlations should modify the normal state behaviour 
between the transition temperature and the mean- 
field transition temperature T*. Early evidence for this 
pre-cursor pairing was the comparison between numeri- 
cal simulations ^ of the attractive Hubbard model with 
the suppression of the Knight shift and nuclear relaxation 
rate [|| between Tc and room temperature in underdoped 
YBCO. The body of experimental evidence for pre-cursor 
pairing in underdoped cuprates has grown considerably, 
for instance, optical conductivity Q and specific heat 
1^ experiments also indicate a loss of weight for normal 
state carriers between Tc and T* . These measurements, 
involve thermodynamics or two-particle response and are 
insensitive to the symmetry of the pairing. The recent 
angle-resolved photo-emission (ARPES) results have 
now confirmed this picture of a loss of single-particle 
spectral weight in the normal state, and in addition, have 
confirmed a symmetry consistent with d2.2_j^2-pairing. 

Staring from the notion that improved understanding 
of the cuprates requires theories of much stronger pairing 
than we are used to, one is naturally lead to generaliza- 
tions of BCS theory. Arguably the simplest candidate 
for such theories is the attractive Hubbard model which 
incorporates increased fiuctuations due to reduced di- 
mensionality with strong (albeit s-wave) pairing and the 
possibility of small pair sizes. Studies of this model have 
generally focussed on either the phase diagram or the un- 
usual normal state. Moreo and Scalapino ||^ followed on 
earlier work by Scalettar et al to estimate the ordering 
temperature by performing finite-size scaling on the pair 
susceptibility. Their Monte Carlo (MC) simulations on 
systems up to size 10 x 10 gave Tktb ~ 0.045 for [/ = — 4 
at quarter-filling (density n = 0.5, energy in units of t). 
Subsequently, Assad et al |jl^ investigated the onset of 
a non-zero superfluid stiffness ps using MC simulations 
on systems of size 8x8. Very recently, Deisz et al |0] 
used the fiuctuation-exchange approximation ||l^ (FI5X) 
to study the superfluid stiffness and the entropy. Most of 
their data was on significantly larger 32 x 32 systems and 
their FEA results agree well with the MC data. However, 



they observe a broad onset of ps (T) (that does not scale 
as expected with system size) instead of the discontinuity 
of the KTB transition. ||l^ The dilemma is that with- 
out a better understanding of the correlation length it 
is difficult to disentangle finite-size effects, and spurious 
long-range order from the expected algebraic order. 

The other class of studies, originated with the MC cal- 
culations of Randeria et al j^] of the temperature de- 
pendence of the spin susceptibility Xspin(T'), the nuclear 
relaxation rate Ti(T) and the density of states at the 
Fermi surface No{T). Their results agreed qualitatively 
with experiments on underdoped cuprates and at an early 
stage suggested that precursor pairing is manifested in 
the normal state of these materials. Further investiga- 
tions in this vein included the MC simulations of Singer 
et al |l^] who calculated the frequency dependence of 
the density of states N{uj). In addition to the poten- 
tially exact MC calculations our earlier gaussian fluctu- 
ation studies ||l|,^ were extended to the self-consistent 
T-matrix approximation (STA) which although ap- 
proximate, allows for the study of different hamiltonians 
and larger systems. From a diagrammatic point of view, 
the STA includes a subset of the diagrams in the FEA. 
Using this self-consistent T-matrix approximation we re- 
cently investigated pre-cursor pairing in a model with 
d-wave symmetry that led to very interesting predictions 
tested in recent ARPES experiments. Q 

Our goal in this paper is to study the critical behaviour 
of the attractive Hubbard model in the parameter range 
relevant for the cuprates, and investigate any connections 
between the Kosterlitz-Thouless-Berezinskii (KTB) tran- 
sition and the anomalous normal state properties. We 
demonstrate that the self-consistent T-matrix approxi- 
mation is sufficiently accurate to produce XY critical 
scaling in the correlation length and pair susceptibility. 
The model has a large critical regime with f (T) falling 
off exponentially from Tktb until ^ ~ 1 (in units of the 
lattice spacing). This happens somewhat before T* and 
above this temperature the pre-cursor pairs are essen- 
tially uncorrelated. Given that the STA includes a subset 
of the FEA diagrams we conclude that both these approx- 
imations capture the delicate interplay p9[ | between the 
carriers, the precursor pairs and the unpairing of vortices. 

The attractive Hubbard model hamiltonian is 



1 



where the first sum runs over pairs of nearest neighbour 
sites on a two-dimensional (2D) square lattice. By vary- 
ing the attraction U one can study the crossover from 
BCS superconductivity in weak coupling, to the Bose 
condensation of local pairs in strong coupling. In 2D, the 
broken U{1) symmetry belongs to the XY universality 
class and one expects a KTB-type transition at non-zero 
temperature. On the other hand, at half-filling, super- 
conductivity becomes degenerate with a charge-density- 
wave instability, and this larger SU (2) symmetry requires 
the ordering temperature to vanish. 

We study this model by self-consistently solving for the 
dressed propagator and four-point vertex. The Green's 
function is defined in terms of the self-energy G{]s., ikn) — 
[ikn — e{k) + fJ, — E(k, ifc„)] . The exact self-energy is 
formally defined through the coupled equations 



I](fc) = UJ2Gi'l-k)Gip)Giq-p)TpA'l) 



(2) 



and 



^k,k' (q) = Ik,k' (q) - ^ Ik,p{q)G{p)G{q-p)Tp^k' (q) (3) 

where the sum (with proper factors) represents a trace 
over intermediate momenta, frequencies and, where ap- 
propriate, spins. Here Ik,k'{q) is the irreducible particle- 
particle vertex, k — (k, ifc„) and q = {(i,iqn) in a tem- 
perature formalism with Fermi and Bose Matsubara fre- 
quencies k„ = (2n + 1)ttT and (/„ = 2mrT respec- 
tively. Determining the full irreducible particle-particle 
vertex is a formidable task and we approximate Ik,k' (q) 
by the bare Hubbard potential. Our description then 
reduces to the self-consistent T-matrix approximation. 
Now Tk^k'iq) — r((7) and (||) simplifies to the separable 
Bcthe-Salpeter equation with solution 



u^njq) 
1 + un{q) 



and (H) reduces to 



S(fc)=5]G(q-fc)r(q). 



(4) 



(5) 



We have subtracted the Hartree shift and the particle- 
particle bubble H(g) = J^k G{q — k)G{k) is defined in 
terms of dressed propagators. 

This amounts to keeping the subset of particle-particle 
diagrams in the fluctuation exchange approximation. 
The self-energy in the FEA includes diagrams in three 
channels 



o 



o 



while we only keep the particle-particle channel repre- 
sented by the first term. Our results confirm that this 
first class of diagrams are relevant and the particle- 
particle contributions are sufficient JTot to produce the 
correct critical behaviour for the densities we investigate. 
At half-filling neither the STA or the FEA methods pre- 
serve the degeneracy between pairing and CDW chan- 
nels that is required to restore the full SU (2) symmetry. 



However, our calculations are sufficiently far from half- 
filling, such that this is not an issue. In fact, our STA 
results agree quite well with the Monte Carlo simulations 
1^ , yielding a suppression of the density of states at the 
Fermi level that tracks the reduction of spin susceptibility 
similar to the Knight shift experiments. A more detailed 
comparison with MC will be presented elsewhere. |Q 

Our goal here is to elucidate the nature of the normal 
state pairing correlations. Specifically we wish to inves- 
tigate the equal-time pair correlation function 



$(r) = (A(0,0)A(r,0)). 



(6) 



This is the Fourier transform of the generalized pair sus- 
ceptibility, which in our approximation has the simple 
form Xpair(g) = n(g)/[1 + C/H(g)]. 

To evaluate Xpair('7)i we employ the standard 
STA/FEA procedure of iteratively solving for (^) and (||) 
on an L X L lattice using fast Fourier transforms. We typ- 
ically keep 256 Matsubara frequencies and to avoid finite- 
size effects study systems of size L = 128. Since this pro- 
cedure uses periodic boundary conditions in space, <I>(r) 
only exhibits bulk character up to distance r ~ L/2. In 
Figure |l| we show our results for <I>(r) at four tempera- 
tures for [/ = — 4 at quarter-filling. We use these param- 
eters for comparison with the MC data |^ which show a 
reduction of the spin susceptibility below T = 0.5. The 
correlation function indeed decays exponentially and at 
large distances shows a very rapid increase upon cooling. 
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FIG. 1. Equal-time pair correlation function <^{r) vs. po- 
sition r — I (a::, 0)1 in units of the lattice spacing for U = —4, 
at quarter-filling, and temperatures T = 0.35, 0.20, 0.15 and 
0.11. 

The correlation length ^ at each temperature is defined 
through the scaling relation |]l3| 



$(r) 



(7) 



at large distances. To extract £,{T), the power law cor- 
rection is small and can be neglected. (For the XY uni- 
versality class the exponent has the theoretical value rj = 
0.25 which indeed is verified below). The temperature- 
dependence of ^ (T) for the same attraction and density 
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is given as the circles in Fig ure The soUd hne is a fit 
to the XY scaUng form [|l3| 

C(T)~exp(A/v/^:^) (8) 

with f3 — I /T the inverse temperature and (3c = I /Tkt b 
the inverse of the KTB transition temperature. The best 
fit yields Tktb = 0.049(2) for = -4 and n = 0.5. For 
comparison we also show the power scaling law forms 
^ {T — Tc)~^ for mean- field and 3D XY exponents 
V = ^ and V — ^ respectively. 
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FIG. 2. Pair correlation length (,{T) vs. temperature T 
for U — —4, at quarter-filling. The solid line is a fit for the 
exponential XY scaling form. The broken and dash lines 
represent power-law scaling for the mean-field and 3D XY 
exponents ~ ^ and |. 

The correlation length clearly scales according the the 
XY critical scaling form and upon heating decreases 
much faster than any power law. It is clearly natural 
to interpret this scaling as due to vortex-pair unbinding 
in normal state. 

Furthermore, the numeric data shows critical scaling 
over a very large critical regime, from Tktb ~ 0.05 up 
to r w 0.30. In fact, we observe critical scaling until 
the correlation length has decayed all the way down to 
a single lattice spacing! The MC simulations and our 
STA results show the suppression of the spin suscepti- 
bility starts to develop below T 0.5. Since ^ does 
not exceed one lattice spacing in the temperature range 
0.35 < T < 0.5 we argue that to describe the onset of 
the pseudogap behaviour one needs to incorporate ampli- 
tude fluctuations in addition to phase fluctuations. The 
size of the critical regime is maximal (in the sense that 
it extends down to ^ = 1 in units of the lattice spacing) 
and the very large scaling regime is reminiscent of the 
quantum critical point. 

A further test for XY scaling is to consider the uni- 
form, static, pair susceptibility which, in our approxima- 
tion, is proportional to T{q) with q = in Eq. ^ Thus 
Xpair diverges as A = —Ull{0) 1. While A reaches 
the value one at a non-zero temperature in mean-fleld 
theory, more accurate calculations such as MC, FEA or 
STA should never give A = 1 for a finite system. In- 



deed, our calculations agree with the FEA calculation of 
Ref. yielding A < 1 for all temperatures. In fact, 
this even holds in three dimensions clearly due to the 
absence of phase transitions in finite systems. However, 
provided we keep the system size L > ^(T), we can ob- 
serve the approach to the divergence in Xpair as Tktb 
is approached from above. In Figure ^ we plot Xpaii (T) 
vs the temperature-dependent correlation length ^(T) to 
the power 2 — rj for U — —4 and n — 0.5 at several 
temperatures. 
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FIG. 3. Pair susceptibility Xpair(r) vs. power- law scaling 
of the pair correlation length ^(r)^"'' for rj = j. 

We have explicitly set rj = j - the universal XY value. 
A fit to the lower temperature data gave rj « 0.248. The 
error in Xpair is larger due to the denominator of (|[) as 
A « 1 at low temperatures. None the less, our data shows 
XY critical scaling over nearly three orders of magnitude. 

The degree to which both ^(T) and Xpair displays 
XY scaling is remarkable. This demonstrates that the 
self-consistent T-matrix approximation which includes 
self-energy corrections but neglects vertex corrections 
still captures |l^] the 'relevant' fiuctuations. Further, 
the very large critical regime in the attractive Hubbard 
Model vindicates our assertion of the importance of small 
distance physics (which is neglected in low-energy, long- 
wavelength expansions). 

We have been maintaining for some time ||l|,^,^ that 
as a model of a 2D short coherence length superconduc- 
tors, the unusual normal state properties of the attrac- 
tive Hubbard Model is of relevance to the high- Tc super- 
conductors. This is especially true for the underdoped 
cuprates which are much more [ p3[ two-dimensional than 
their optimally-doped or overdoped counterparts. 

It is therefore natural to conclude that our results have 
far-reaching implications for the cuprates. The picture 
that emerges is that upon cooling the normal state be- 
haviour is strongly infiuenced by the 2D XY critical point 
over a large temperature range, until three-dimensional 
coherence starts to develop just above the observed (3D) 
critical temperature. Thus, the strong 2D pairing corre- 
lations imply that (zero-field [|^) vortex-pair unbinding 
is important in the normal state of the cuprates and that 
vortices provide an additional (strong) scattering mech- 
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anism. This scattering channel disappears below and 
although this s-wave model of course does not contain 
the d-wave pairing observed in ARPES experiments, this 
framework is clearly consistent with the observation of 
broad spectral peaks in the normal state which sharpen 
dramatically below Tc. The short distance over which 
phase-coherence is established then also explain why very 
strong paraconductivity is not observed over a large tem- 
perature range. 

These results also emphasize that in order to study the 
onset of order numerically, one needs very large systems 
to eliminate finite size effects. The transition observed by 
Deisz et al is therefore indeed a KTB type transition 
although the 32 x 32 systems they focussed on may be 
on the borderline of being influenced by finite-size effects. 
Whereas the STA/FEA methods may have problems near 
half- filling where the model has SU{2) symmetry, these 
methods are very successful at keeping track of relevant 
pairing correlations away from half-filling while ensuring 
that large correlation lengths do not mimic spurious long- 
range order. These methods should therefore be used 
in studying a variety of models to complement the MC 
calculations which are restricted to smaller lattices, but 
of course have other advantages. 

It should be noted that several other groups have em- 
phasized 2D pair correlations in the cuprates, dating 
back to Uemura's |2^] investigation of the superfluid den- 
sity. Also, the physics we discuss here and introduced in 
Refs. [p],p p^jl7| ] bears some resemblance to the frame- 
work developed later by Emery and Kivelson |p6|| - al- 
though they emphasize the proximity of an insulating 
state in the cuprates. The attractive Hubbard model is, 
however, also susceptible to a metal-insulator-transition 
near half-filling due to the charge-density-wave correla- 
tions, although this simple model may not capture the 
doped Mott insulator physics of the cuprates. 

Above we have demonstrated the connection between 
XY fluctuations and pseudogap physics. It is important 
to emphasize that it is crucial to also be in the intermedi- 
ate coupling (short coherence length) regime of the phase 
diagram. In the weak coupling limit XY fluctuations do 
not produce pseudogap behaviour over an appreciable 
temperature range. [ P7[ 

To summarize we have demonstrated that the self- 
consistent T-matrix approximation can capture the nor- 
mal state XY physics of the unbinding of vortex pairs 
for the attractive Hubbard model in 2D. The scaling 
regime is very large, extending beyond 5 times Tktb, un- 
til the correlation length has decreased down to a single 
lattice spacing. The large regime influenced by vortex- 
pair unbinding suggests that the excitations may be rel- 
evant to the normal state pseudogap behaviour of under- 
doped cuprates. The large window of temperatures where 
^(T) ^ 10 also implies that when investigating the devel- 
opment of order, it is important to complement MC cal- 
culations with STA/FEA calculations on large systems. 
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